Gravitational Radiation from Primordial Helical Inverse Cascade MHD Turbulence 
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We consider the generation of gravitational waves by primordial helical inverse cascade magneto- 
hydrodynamic (MHD) turbulence produced by bubble collisions at the electroweak phase transition. 

■ We extend the previous study [l|] by considering both currently discussed models of MHD turbu- 
' lence. For popular electroweak phase transition parameter values, the generated gravitational wave 
, spectrum is only weakly dependent on the MHD turbulence model. Compared to the unmagne- 

■ tized electroweak phase transition case, the spectrum of MHD-turbulence-generated gravitational 
5^ ' waves peaks at lower frequency with larger amplitude and can be detected by the proposed Laser 
CL(' Interferometer Space Antenna. 

\o ■ 
<N : 

I ^- I. INTRODUCTION 

^ : 

The direct detection of a gravitational wave background will provide us with a probe of physical conditions in the 
Q . early Universe at the epoch when the gravitational radiation was generated. This is because gravitational waves are 
^ ■ massless and after generation weakly-coupled gravitational radiation propagates freely. Hence, once generated, any 
c/3 , gravitational wave spectrum retains its shape, with all wavelengths simply scaling with the expansion of the Universe. 
, ^ , • Gravitational wave astronomy holds the potential for helping construct a picture of the Universe at energy scales even 

higher than those the Large Hadron Collider will reach, to probe in detail physics at the electroweak energy scale. 
CO ' There arc various mechanisms that might generate gravitational waves in the early Universe, for reviews see 0,0]. 
^ ' A well-known one is parametric amplification of quantum fluctuations Q, which can also take place shortly after 
0^ inflation Q . Other mechanisms include bubble wall motion and collisions during phase transitions; gravitational wave 
generation during the electroweak and QCD phase transitions are discussed in Refs. [ 9 | . [lOll . Cosmic strings and 

other defects can produce relic gravitational waves [ll| ■ Cosmological magnetic fields and hydrodynamical 

or magnetohydrodynamical turbulence can also induce primordial gravitational wavesfll. 1 1 5l. 1 1 6l. [itI . [18| . 

Gravitational waves will also provide information on symmetries present in the Universe. In particular, the relic 
gravitational wave background carries information about parity asymmetry, that might be generated at the electroweak 
' phase transition, or earlier. Parity symmetry testing based on the direct detection of gravitational waves in the near 
future (with currently planned missions) seems to be promising for gravitational waves with frequencies around the 
Laser Interferometer Space Antenna (LISA) sensitivity band at 0.1 — 100 mHz [l9|. The Standard Model electroweak 
phase transition, when parity symmetry breaking and magnetic helicity production might be expected [20j, produces 
^ gravitational waves with a significant amplitude in this frequency range. Recently, some of us have investigated the 
5^ ', generation of gravitational waves through inverse cascade MHD turbulence [l| and found that even a small amount of 
initial helicity enlarges the range of phase transition parameters for which gravitational waves are potentially detectable 
by LISA. This gravitational wave signal has an important feature: it is circularly polarized. For direct cascade 
hydrodynamical turbulence some of us have previously estimated the polarization degree (l7j . which is obviously 
turbulence model dependent. Polarized gravitational waves are present in other models [2]| . and the polarization of 
the gravitational wave background is in principle observable, either directly [l^l or through the cosmic microwave 
background ^,|2l|. 

The detection prospects of gravitational radiation from the early Universe depend on the energy scale when gravi- 
tational waves were generated and even more so on the source duration time and efficiency of the generation process 
itself. In particular, to produce a detectable signal the electroweak phase transition must be strong enough HI, HI] ■ 
The currently popular standard model of particle physics does not have a strong enough electroweak phase transition. 
Also, all currently discussed electroweak phase transition models do not predict an observable gravitational radiation 
signal if one assumes only bubble collisions [26| . It is of particular interest to study gravitational wave production 
duringphase transitions in the minimal and next to minimal supersymmetric standard models (MSSM and nMSSM 
[13)i H- The gravitational wave signal produced in the MSSM model is below the LISA sensitivity, while the nMSSM 
results in a stronger electroweak phase transition 27] and as a consequence a detectable gravitational wave signal Q. 
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On the other hand, sources associated with stochastic vector fields in plasma, i.e. kinetic velocity or magnetic 
fields, significantly change gravitational wave detection prospects [l^. In the case of unmagnetized hydrodynamical 
turbulence the peak frequency of the gravitational wave power spectrum is determined by the characteristic time of 
turbulence, i.e. the inverse turn-over time of the largest eddy. Another important characteristic is the energy scale 
when gravitational radiation is generated. In recent modifications of the standard electroweak phase transition model 
where the transition is moved to a higher energy scale [25l [29l | , the gravitational wave power spectrum peak frequency 
is shifted to a higher frequency, which, since the gravitational wave spectrum is sharply peaked, reduces the possibility 
of detection by LISA. The presence of a cosmological magnetic field affects the dynamics of turbulence. In the case 
of MHD turbulence the presence of an energy inverse-cascade leads to an increase in the effective size of the largest 
eddy (now associated with a helical magnetic field), and can result in the gravitational wave power spectrum peaking 
in the LISA band, with amplitude large enough to be detected by LISA 1]. 

In this paper we extend the study of Ref. [I| by also considering a different model of MHD turbulence. In particular, 
Ref. [H adopts the inverse cascade model of Refs. [13, HH, and in this paper we also include in the consideration the 
model of Refs. [H, IH, Is^]. In both cases, gravitational wave generation is considered assuming non-zero magnetic 
helicity during the phase transition. 

For our analysis in this paper we adapt the technique developed in Ref. [3l ■ We model MHD turbulence and obtain 
the gravitational wave spectrum by using an analogy with the theory of sound wave production by hydrodynamical 
turbulence [ssl. [36l. [stI. [Sq . We perform the computation of the gravitational wave spectrum in real space, instead of 
using conventional Fourier space techniques as in Refs. [^.[lH. This makes the physical interpretation of all quantities 
straightforward. We use natural units with h ~ c = ks = I throughout. 

The paper is organized as follows. In Sec. II we summarize the general formalism of gravitational wave generation. 
In Sec. Ill we describe the turbulence models we consider. In Sec. IV we derive the gravitational wave spectra, and 
we conclude in Sec. V. 



II. GRAVITATIONAL WAVE GENERATION FORMALISM 



We assume that primordial turbulence is generated at time t^, — corresponding to temperature at the phase 
transition — at characteristic proper length-scale = 27: /k^ (where k^ is the corresponding proper wavenumber) 
and with characteristic velocity perturbation vq. We assume that the kinetic and magnetic Reynolds numbers are 
much greater than unity on scales of order ^q, otherwise there is no turbulence. We also assume that the duration 
of the turbulence, tt, is no larger than the Hubble time at the phase transition H~^, where H^, is the value of the 
Hubble parameter at the pha se transition, so the expansion of the Universe may be neglected during the generation 

16 . This adiabatic assumption will be valid for any turbulence produced in a realistic 
3£ |. Therefore, the gravitational radiation equation in real space can be written as (40| 



of gravitational radiation 



cosmological phase transition 



g2 

V2/iy(x,<) - — /iy(x,i) = -16^G5y(x,t), (1) 



where /iy(x, i) is the tensor metric perturbation, t is physical time, i and j are spatial indices (repeated indices 
are summed), and G is the gravitational constant. Sij is the traceless part of the stress-energy tensor Tij{x,t), 
(constructed from either kinetic or magnetic turbulence vector fields), given by (40| 

5.,(x,i) = T,,(x,i)-i%r,^(x,i). (2) 

Since the turbulent fluctuations are stochastic, so are the generated gravitational waves. Our goal is to derive the 
energy density spectrum of these gravitational wave perturbations at the end of the turbulent phase; after that the 
amplitude and wavelength of the gravitational radiation scales simply with the expansion of the Universe. The energy 
density of gravitational waves is [3[ 

- ^(a.%(x,t)a./.,(x,^)) . -y d3x-dV- ^ |x-x'||x-x1 ' 

where the brackets denote an ensemble average over realization of the stochastic source, and the times = t — |x — x'|, 
t" = t- |x-x"|. 

We consider metric perturbations in the far-field limit (i.e. for x ^ d, where d is a characteristic length scale of 
the source region), where gravitational waves are the only metric perturbations [40| . and replace |x — x'| by |x| in 
Eq. ([3]). For the stationary turbulence we consider here, we define the gravitational wave spectral energy density. 
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/(x, Lu), as the one-dimensional temporal Fourier transform of the autocorrelation function of the temporal derivative 
of the tensor metric perturbations i(x, t), 

/(x,c.) = i-ydre'--L(x,r), (4) 

where uj is the angular frequency and 

L{x,t) = — i— {dth^j{x,t)dth^j{yi,t + T)). (5) 
Accounting for Eq. ([3]) , the gravitational wave energy density is 

PGw(x) y"dw/(x,w). (6) 



Following Ref. 'l8j it can be shown that 



47r^w^G'w^ f / , X 



/(x, iu) = — / d^x'i?,,,, x', — a>, , (7) 



X 



where k is a proper wavevector and w = p + p, with w, p, and p being respectively the enthalpy, pressure, and energy 
density of the plasma. The quantity Hijijix.' ,]<.,ll!) is the double trace of the four-dimensional temporal and spatial 
Fourier transform^, taken with respect to x" — x' and r, of the two-point time-delayed fourth-order correlation tensor 



That is, 



if,,fci(x',k,c.) = J d^l dr e'(-^-'^ «)i?,,feKxM,r), 



(9) 



where ^ = x" — x' . 

Given the function iJyfc; (x', k, w), we may use Eqs. ([6]) and ([7|) to compute the gravitational wave energy density, 
assuming that the source is statistically homogeneous, so that the averaged correlators of the stress-energy tensor 
have no spatial dependence, and isotropic, so that the correlation between two spatial points depends only on the 
distance between the points and not on the direction. With these assumptions, Eq. ([9]) reduces to 

H,,M{^',k,Lu) = H,,m{Klo) = ^ j A^i dr e'^^^-'^-^^ R,,ki{i,T), (10) 

so Hijki{'x.w,w) — Hijki{uj,uj) is independent of the observation direction x, as expected on physical grounds. 

Finally, since the gravitational waves propagate freely after been generated, the expansion of the Universe is 
accounted for by a simple re-scaling of the frequency and the amplitude by a factor equal to 

«*~.8xlo-fH^vH2V'^ (11) 



ao \ Ti, J \ g- 

where 5^ is the number of relativistic degrees of freedom at the temperature T^,. For the standard model degrees of 
freedom we have 5* = 106.75 as T — > 00. With these assumptions and considering a stochastic turbulence source 
lasting for a finite duration tt, the total gravitational radiation energy density spectrum at a spatial point and a 
time can be obtained by integrating over all sources within a spherical shell centered at that observer, with a shell 
thickness corresponding to the duration of the phase transition, and a radius equal to the proper distance along any 
light-like path from the observer to the source. Then we get [ll, [l^ 



PGiv(w) = '^P^'^ ^ 16T:^uj^Gw'^TTH^jij{io,io), (12) 

CtlHUJ 



^ Consistent with current observational indications (4ll |. we assume flat spatial hypersurfaces and use spatial Fourier transform relations 
Fj(k) = f (fxe^^ '^Fjix) and Fj{x) = J ifk e-'''-^Fj (k)/(27r)3 . 
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where w is the angular frequency measured at the moment of generation of the gravitational radiation. As shown 
in Ref. p^ . while Hijij{uj,uj) is a complicated function of w, the resulting gravitational wave signal strength can be 
simply estimated with ±25% accuracy by working in the k ^ limit and replacing Hijij (k, to) by the corresponding 
Hijij{Q,oj) — this is the main point of the aero-acoustic approximation initially used for sound wave generation by 
turbulence [361] . It can be shown that this approximation works rather well for low velocity turbulence, with Mach 
number M less than unity, i.e. for M — vq < 1, where vq corresponds to the turbulence velocity of the largest eddy. 

We assume that turbulence is caused by bubble collisions during the phase transition. In general, phase transitions 
can be described by two key parameters: a = Pvac/pthcrmai, the ratio of the vacuum energy density associated with 
the phase transition to the thermal energy density of the Universe at the time (this characterizes the strength of the 
phase transition); and /3~^, which sets the characteristic time scale for the phase transition, as well as the temperature 
T-i, when the phase transition occurs. All other parameters, such as the efficiency factor k which gives the fraction of 
the available vacuum energy which goes into the kinetic energy of the expanding bubble walls (as opposed to thermal 
energy), the bubble wall velocity wt,, and the turbulent eddy velocity vq, can be expressed through T*, a, and 
In particular 0, [13, SI 



vt{a) = 
vo{a) = 



1 + 0.715a 



0.715a 



4 

27 



3a 
T 



1/2- 



l/V3+(a2+2a/3)i/2 



3Ka 
+ 3Ka 



1 + a 

1/2 



(13) 
(14) 
(15) 



The parameter is determined from the bubble nucleation rate per unit volume F = Foe'^* [3] where Fq is determined 
by initial conditions. For the standard model = lOOi?^. The temperature T,, and the Hubble parameter at the 
phase transition H^, are related through 



2 _ STr'^GTt 
90 ■ 



(16) 



In the next section we discuss two specific models for the turbulent motions. 



III. TURBULENCE MODELS 



In this paper we are interested in the generation of gravitational waves by helical MHD turbulence, that is turbulence 
in a magnetized plasma with non- vanishing magnetic helicity, Ti.M, defined by 



HAiit) = J d^x A-B, 



(17) 



where A is the vector potential of the magnetic field B. Magnetic helicity is odd under discrete P and CP trans- 
formations, so the presence of magnetic helicity in our Universe would be an indication of macroscopic P and CP 
violation. 

While MHD turbulence is isotropic on large scales, it is locally anisotropic on small scales [4^, resulting in small- 
scale anisotropy in the generated GW background. However, GWs are generated mainly by the largest eddies llSl so 
we adopt an isotropic turbulence model with velocity field and magnetic field two-point correlation functions |35| 



{vnk,t)v,ik',t + r)) = 6ik - k') Fr^{k,t) f[v{k),T] 



(18) 



(&*(k, t)b,{k', t + r))= 5{k - k') Flf(k, t) f[v{k),T] 



(19) 



Here v is the turbulent- motion velocity field, b = B/V47rw is the characteristic magnetic field perturbation B 
expressed in velocity units, and 



F;;.(k,r)=Py(k) 



Ey{k,t) 

47rfc2 ' 



(20) 



(21) 
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In the above equations the projection operator -Pij(k) = 5ij — kikj/k^, 6ij is the Kronecker delta, k = |k| 
Eiji is the totally antisymmetric tensor, and r](k) is an autocorrelation function that determines the character- 



istic function / 

f[r]{k),T] = exp 



r]{k),T] describing the temporal decorrelation of turbulent fluctuations. In the following we use 
7r7y^(fc)r^/4] [4J|. The function Ey{k,t) in Eq. (|20p is the so-called kinetic power spectrum and is 



related to the kinetic energy of turbulence through 

Ey{t)^ JdkE,{k,t), (22) 

while in Eq. (j2ip EMik,t) and HM{k,t), the magnetic field energy and magnetic helicity power spectra, denote the 
symmetric and antisymmetric parts of the two-point magnetic field correlator. They are related to the magnetic 
energy and magnetic helicity densities through 



£M{t) = I dkEM{k,t) (23) 

and 

Huit) - Idk HM{k,t), (24) 



respectively. 

It is worth noting that for all magnetic field configurations, the magnetic helicity spectrum must satisfy the "realiz- 
ability condition" [4^: |iJA/(fc,t)| < 2_Ej\/(fc, t)/fc. Introducing the magnetic correlation length (i.e., the characteristic 
proper length associated with the large magnetic energy eddies of turbulence), 

Jdkk~^EM{k,t) 

Ul(t) = TT-TTx , (^5) 

CM[t) 

the integral form of the realizability condition is 1 7^1/(0 1 ^ '^£.M{t)£M{t)- In the following, we assume that a magnetic 
field with (positive) magnetic helicity has been generated during a phase transition through some mechanism [20| . 
We also assume "small" initial magnetic helicity, in the sense that the parameter 

^ Ti.M{t*) , . 

^* 2S,m{U)Sm{U) ^ ' 

is much smaller then unity. 

After generation, primordial helical turbulence decays freely. While it is widely accepted in the literature that the 
free decay of helical MHD turbulence is a two-stage process, the details of the two stages are still under discussion. 
Roughly speaking, in the first stage, the system proceeds through selective decay of magnetic modes: magnetic 
power on small scales is washed out by turbulence effects more effectively than on large scales. During this process 
the magnetic correlation length grows while the magnetic energy decays in time. The first stage ends when quasi- 
conservation of magnetic helicity starts to trigger an inverse cascade of the magnetic field: small-scale modes are 
no longer completely dissipated by turbulence, rather part of their energy is now transferred to larger scales. This 
causes a faster growth of the correlation length and a slower dissipation of the magnetic energy in comparison to the 
non-helical case. 

Different analyses in the literature seem to agree on the details of the first stage of direct cascade. The details of 
the inverse-cascade regime, however, are still being debated and there are two main models, that of Refs. [sol l3l|. 
here referred to as Model A, and that of Refs. i33i,i34,], which we will call Model B. In the following, we consider both 
models. 

A. First stage: direct cascade 

In both models the dynamics of the first decay stage is governed by a direct cascade of magnetic energy density 
lasting for a time Tsq = sqTq, a few times (sq ~ 3 — 5) longer than the characteristic largest-eddy turn-over time 
To = ^o/vq — 27r/fcofo- ^ At the end of this stage equipartition between kinetic and magnetic energies is reached [45}. 



^ To be more precise, from the results of Refs. [33l. l34|. the time when the system enters the inverse-cascade regime (in Model B) could 
depend on the fraction of the maximal initial magnetic helicity f*. However, those results were obtained assuming an initial magnetic 
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During the first stage magnetic energy density flows from large to small scales and finally dissipates on the scale 
Id = '^T^/kd (kd ^ ko) where one of the Reynolds numbers becomes of order unity. Due to the selective decay 
effect [43 magnetic helicity is nearly conserved during this stage [sil [s^l- To compute the gravitational waves 
generated by decaying MHD turbulence, we assume that decaying turbulence lasting for time Tso is equivalent to 
stationary turbulence lastin g fo r time Tso/2. This can be justified using the Proudman argument for (unmagnetized) 
hydrodynamical turbulence |35l. Issj] . Consequently, when computing the induced gravitational waves we ignore the 
time dependence of EM{k,t) and HM{k,t). Also, for EM{k,t) and ri{k) we use the Kolmogorov model (ioj with 

EM{k,t) = CKe"''k-'''\ (27) 

and 

7;(fc) = ^ (28) 

both functions being defined over the range of wavenumbers k^ < k < kd- Here Ck is a constant of order unity and 
e ~ fco^^o ^he energy dissipation rate per unit enthalpy. Taking into account that equipartition between kinetic 
and magnetic energy densities is established and maintained during Kolmogorov turbulence, the total (kinetic plus 
magnetic) energy density of turbulence is simply double the expression in Eq. (|27p . 

We note that while turbulence in the early Universe has to be relativistic, we have assumed a Kolmogorov spectrum 
for both b and v, which is valid for non-relativistic turbulence. Based on arguments of Ref. p^ . however, we expect 
that our estimate gives the correct qualitative features of the resulting gravitational radiation spectrum. 

B. Second stage: inverse cascade 

1. Model A 



Model A has been discussed in Ref. Jl|. Here, we summarize those results. At the end of the first stage turbulence 

relaxes to a maximally helical state [3ll . [33j . Accounting for conservation of magnetic helicity, the characteristic 

1/2 1/2 

velocity and magnetic field perturbations at this stage are ui ~ C^' vi^ and bi ~ C,J 60 ■ The dynamics of the second 
stage is governed by a magnetic helicity inverse cascade. If both Reynolds numbers are large at the end of the first 
stage, magnetic helicity is conserved during the second stage. The magnetic eddy correlation length evolves as [sol-lsij 



CA/(t)~;o(^l + -j , (29) 

1 /2 

where ri ~ Iq/vi = to/(J is the characteristic eddy turn-over time at the beginning of the second stage. The 
magnetic and kinetic energy densities evolve as (3Q . i3lj 



SMit) ~ w6^( 1 + - 

Tl 



-1/2 



-1 



£^{t) ~ wvf(l + —] . (30) 



"^1 . 

These imply that the characteristic turn-over (rto) and cascade (rcas) timescales evolve as 

Tto ^ Teas ~ n ^1 + -^^ . (31) 

To compute the gravitational waves emitted during the second stage we use the stationary turbulence model that 
has the same gravitational wave output. Introducing the characteristic proper wavenumber k^{t) — 271/^^1(1) and 



spectrum of the form BmC^iO) oc with p a positive real number, while in our case the initial spectrum is a red spectrum of the 
Kolmogorov type. In the absence of firm results on this matter, and in order to avoid inessential complications, here we assume that in 
both models of turbulence Tsq = sotq with sq ~ few. 
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using Eq. (|30p we find £m — wt;Jfc^(t)/fco since hi ^ vi. Since the kinetic energy density is dissipated more efficiently 
than the magnetic one, we can neglect its contribution in the following discussion. The time when turbulence is 
present on scale ^m(0 is determined by Eq. ([3T|l . which can be rewritten as 



fco 

hit) 



(32) 



So instead of considering decaying turbulence, we consider stationary turbulence with a scale-dependent duration 
time (time during which the magnetic energy is present on that scale), Tgi ~ ri[fco/A:]^ (for k = this coincides with 

"^cas ) ■ 

The expression for £m yields the time-independent magnetic field energy and magnetic helicity power spectra 

EM{k,t) = ^ = i kHMik,t), ks<k< ko. (33) 

/Co Z 

Here Ci is a constant of order unity, ks is the smallest wavenumber where the inverse cascade stops, and the second 
equation follows from saturating the causality condition. For the second stage autocorrelation function, which is 
inversely proportional to the turn-over time ()31|) . we assume 

to be as close as possible to the first stage description where ?7(fco) = ^/^/to- 

At the largest scales there is no efficient dissipation mechanism, so the inverse cascade will be stopped at the scale 
ls{t) = 27r/fc5 where either the cascade timescale Teas reaches the expansion timescale H^^ = H^^{ti,), or when 
the characteristic length scale £,Mi't) — h reaches the Hubble radius. These conditions are ^^^Is/vqIq < H^"^ or 
Is < H^^ (the cascade time is scale dependent and maximal &% k = ks). Defining 7 = lo/H~^ ~ VbH^,(3^^, it is 
easy to see (taking into account that vq, 7, and are less than unity) that the first condition is fulfilled first and 
consequently 

^ < ( ^)'''cy^ (35) 

ks V 7 / 

1/2 

To have an inverse cascade requires ko/ks > 1, leading to a constraint on initial helicity: 7 < MQJ . 



2. Model B 

Instead of assuming the MHD turbulence model of Refs. [s^, HH, in this subsection we consider the alternative 
MHD turbulence model of Refs. [H, [ll] (also see Ref . [l^l ) with different evolution laws for the magnetic correlation 
length and magnetic and kinetic energy densities, 

a/(t)^;of l + , (36) 

8M{t)~whl{l + ^^ ' , (37) 
t 



8,{t)~wvi\l + -j , (38) 

respectively. Taking into account the above equations, it is easy to verify, however, that the characteristic turnover 
and cascade timescales evolve as in Model A. 



^ Taking into account that the turbulent eddies are formed through bubble expansion and bubble collisions, the energy-containing scale 
lo is determined by the bubble wall velocity vi, and the expansion time as Iq ~ vi,[3~^. We note that 7 determines the number 
(A^eddy — 7~^) of turbulent eddies within the Hubble radius H^^ . 
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It should be noted that, although the magnetic and kinetic energies scale in time in the same way, equipartition is 
not generally reached [ssj . However, if the initial ratio of kinetic and magnetic energies is taken to be of order unity, 
equipartition is established and lasts for all time. Applying these results to the case at hand, after the first stage of 
direct cascade, at the end of which equipartition between magnetic and kinetic energy is established, these follows a 
stage of inverse cascade which preserves equipartition of energy. Hence we conclude that kinetic energy cannot be 
neglected in the second stage for this particular model of turbulence. In fact the amount of gravitational radiation 
produced in the second stage (which, as we will see, is much greater than that generated in the first stage) is exactly 
twice that produced considering only magnetic energy. 

Proceeding as for the discussion for Model A between Eqs. (|3ip and (l33l) . we find that the relation ~ 
wvf [A:^(<)/fco] does not change, whereas for the cascade timescale we get 

-| 3/2 

(39) 



The expression for the time-independent magnetic field energy spectrum does not change, however, now it is valid for 
ks < k < ko with ks depending on the particular model of turbulence adopted (see below). 
In Model B, the autocorrelation function is 



Tl \ko 

while the conditions that determine when the inverse cascade stops are CvT^^^'l^^/wo^o^^ < or Is < H^^. However 
in this case it is easy to see that the former condition is fulfilled before the latter one, and, consequently, we get 



ks V 7 

The condition to have an inverse cascade, kg/ks > 1, leads to the same constraint on initial helicity found in Model 
A, i.e. 7 < A/C^ ■ 



C. Stress-energy tensor and source for gravitational waves 



The magnetic field perturbation stress-energy tensor is 

(x,t) =w&,(x,t)6,(x,i). (42) 

For the first decay stage we compute for this magnetic part and then double the result to account for approximate 
magnetic and kinetic energy equipartition for Alfven waves. 

To compute the function Hijij (k, ui) , we assume Millionshchikov quasi-normality (35| and adopt the (k — s- 0) aero- 
acoustic approximation, which is accurate for low Mach numbers and slightly overestimates the gravitational wave 
amplitude for Mach numbers approaching unity (for details see Sec. HI of Ref. [3). The final result is [isj 

Hm (k, ^) ^ H,,,, (0, uj) = J^^ — exp i^-^j^j erfc(--^y^ j . (43) 

Here, erfc(x) is the complementary error function defined as erfc(a;) = 1 — crf(a;) , where erf (x) — J^^ dy exp(— y^) is the 
error function 47]. The integral in Eq. is dominated by the large scale (fc ~ fco) contribution so, for direct-cascade 
turbulence during the first stage, the peak frequency is [l8| 

u;\il 0, koM. (44) 
To compute the gravitational wave source during the second stage we consider Model A and Model B separately. 



1. Model A 



During the second stage, according to Eq. ((30|) . kinetic energy can be neglected compared to magnetic energy. 
Proceeding as in the case of Kolmogorov turbulence and making use of the aero-acoustic approximation, we find 

i?,,,(k,.) ^ F,,(o,.;) - ,,^3/2fc„ ^ F"^rc:^j''''("d^j' ^ ^ 
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where ks can be found by saturating Eq. (p5|) . In this case the integral is dominated by the large scale {k ~ ks) 
contribution and is maximal at 



'^max — ; — ^T"-"* 



2. Model B 



(46) 



In Model B of freely decaying MHD turbulence, as discussed in Sec. II. B. 2, during the second stage kinetic energy 
is approximatively in equipartition with magnetic energy. To compute Hijij (k, uj) , we proceed as for the case just 
discussed for Model A. The only differences reside in the expressions for the autocorrelation function and the smallest 
wavenumber where the inverse cascade stops, kg [whose value can be found by saturating Eq. (I4ip ]. For the case at 
hand we find 

(k, ^) ^ H,j,j (0, u) = , / — exp ( - ■^mjZH ) ^rfc | - I • (47) 



Again, the integral is dominated by the large scale (fc ~ kg) contribution and is maximal at 

-'Ail - " - 2.if. . (48) 

The fact that the frequency at the peak of the generated gravitational wave spectrum is independent of the particular 
model of turbulence [compare Eqs. and (|48p ] is easily understood if one observes that it is the inverse of the 
maximum cascade timescale r^g'^, which is t^'^ — H^^ by definition. 

IV. GRAVITATIONAL WAVE SPECTRA 

The total gravitational wave energy density spectrum pGwl^^) at a given space-time event is obtained by integrating 
over all source regions with a light-like separation from that event, and includes contributions from gravitational wave 
generated during the first and second stages. For the first stage (with duration time t!^^ = sqTq) Pq]^(w) is given by 
Eqs. (21) and (A3) of Ref. 18]. For the second stage contribution we must account for the scale dependence of the 
cascade time. The total GW fractional energy density parameter at the moment of emission is [l^ 

17gw(-.) ^ 105^ ^4")i/,g(0,c..). (49) 

^' m 

Here the index m runs over / and // for the first and second decay stages, cj* = u}{t^) is the angular frequency of the 
gravitational wave at the moment of its emission, and Hq is the current value of the Hubble parameter. 

The current gravitational wave amplitude is related to the current fractional energy density parameter through 

hcU) = 1.26 X 10-i« [hi r!Gw(/)] , (50) 

where /iq is the current Hubble parameter in units of 100 kmsec^^Mpc^^ ,3]. Using the expressions found for the 
Hijij tensor, we obtain 



7, (f\ o in-i4 MOOGeVWlOO\-'^^-- r („) „ATj{m),^ .V'^ 

/ic(/) - 2 X 10 ( — — — 1 ( — 1 Z^[4 ^*-ff*^W(0'^*)J ■ (51) 

Here, / = {ai,/ao)fi, is the linear frequency, with f^, = a;*/27r. 

Figure 1 shows the gravitational wave amplitudes, hc{f), from pure hydrodynamical turbulence (no inverse cascade) 
and for Models A and B for two different values of initial magnetic helicity, = 0.15 (left panel) and = 0.05 (right 
panel). Accounting for inverse-cascade MHD turbulence, the gravitational wave spectrum has two peaks: the first 
higher-frequency one is associated with direct-cascade hydrodynamical turbulence while the second lower-frequency 
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FIG. 1: The spectrum of the gravitational wave strain amplitude, hc{f), as a function of the frequency / for a first-order 
phase transition with — 100, T, — 100 GeV, a — 0.5, and /? = 100//*, from hydrodynamic Kolmogorov turbulence with 
zero magnetic helicity (solid lines) and for the two MHD turbulence models, Model A (dash-dotted lines) and Model B (dashed 
lines). The left panel corresponds to initial magnetic helicity = 0.15, while C* ~ 0.05 in the right panel. In both panels 
the bold solid line corresponds to the 1-year, 5a LISA design sensitivity curve ]49\ including confusion noise from white dwarf 
binaries [4911. 
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FIG. 2: As in Fig. 1, except now T* = 250 GeV. 



one is induced by inverse-cascade MHD turbulence. The arnplitude of gravitational waves emitted during direct- 
cascade unmagnetized turbulence peaks at current frequency [13, HI, [M, 2^ 



HI) ^Mkoa^^(M\ 

./max r,_ „ \ „. I ^ fj ] Jn- 



2tt oq 



Vb 



Here we have used 2Tr/ko = Iq — Vb/S ^, where Vb is the bubble wall velocity, and defined 



1.6 X 10"'' Hz 



ao 



100/ 



100 GeV 



(52) 



(53) 



The second peak is at a lower frequency compared to the unmagnetized case and is independent of the adopted 
turbulence model (since it is determined by the Hubble frequency at the moment of gravitational wave generation). 
In fact, the second peak frequency is 



(11) 

max 



fH- 



(54) 



Finally, as it is straightforward to establish, the amplitude of MHD-turbulence-generated gravitational waves at the 
peak in Model A is about a factor {M/^)^/^C.^^^ larger than that in the unmagnetized case (and thus strongly depends 
on initial magnetic helicity), while it is about a factor \^ {M smaller than that in Model B. 

In Fig. 2 we show the strain amplitude hc{f) for a higher phase transition temperature, T* ~ 250 GeV, with the 
rest of the parameters taking the same values as in Fig. 1. Comparing Figs. 1 and 2 we see that if the phase transition 
occurs at higher energy the peak frequency is higher, closer to the LISA design sensitivity peak. 
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FIG. 3: The LISA sensitivity region for Model A in the P/H^ and T« parameter plane for a phase transition with vacuum 
energy a = 0.1 (left panel) and a = 0.5 (right panel) [2^. The regions for = and = 0.15 coincide at these temperatures 
for Of = 0.1 (left panel). A point in parameter space is considered detectable if at any frequency its value of hc{f) is detectable 
at a signal-to-noise ratio of 5 in a one-year integration, including the confusion noise from white dwarf binaries, based on Refs. 




T, (GeV) T, (GeV) 

FIG. 4: The LISA sensitivity region for Model A in the a and T* parameter plane (left panel) with = 100 and = 100, 

and in the <?* and parameter plane (right panel) with 13/H-^ = 100 and a = 0.1. A point in parameter space is considered 
detectable if at any frequency its value of hc{f) is detectable at a signal-to-noise ratio of 5 in a one-year integration, including 
the confusion noise from white dwarf binaries, based on Refs. [4^ . 



The gravitational wave signal detection prospects depends crucially on several parameters: the temperature (the 
energy scale) at which the phase transition takes place Tt, the number of relativistic degree of freedom g^,, the phase 
transition model parameters a and [3/ Hi,, and the initial magnetic helicity. 

The dependence on Qi, is quite weak, see Eq. (j5ip . It affects both peak frequencies (associated with hydrodynamical 
and MHD turbulence). The initial magnetic helicity plays a more substantial role if the phase transition takes places 
at a higher temperature [1^. The phase transition model determines the value of the parameter a. The tensor Hijij 

appearing in Eq. ([5T|) depends on a as well as on the duration time of the direct cascade stage r-^/'' [ll| . In addition, 
the frequency of the first peak (~ k^M) depends on a through the bubble wall velocity [vb) and other parameters, 
see Eqs. (fT3|) — (|15p . Increasing a makes the gravitational wave signal stronger, see Fig. 2 of Ref. [1^. In the context 
of the detection of gravitational waves, the phase transition energy scale has two different effects: the peak frequency 
is proportional to T^, but the peak amplitude is inversely proportional to T^. These effects must be accounted for 
together, see Fig. 3 of Ref. [1^. Another important parameter is 13/ Hi, which is model dependent and affects the first 
peak frequency, see Eq. (15^. as well as the amplitude of the gravitational wave signal. The dependence on (3 /Hi, can 
be seen in Fig. 5 of Ref. [28|: a larger [3 /Hi, results in a lower amplitude signal with a higher first-peak frequency, 
while the position of the second peak, associated with inverse-cascade turbulence, is not affected. Figure 3 shows the 
LISA sensitivity region in the (3 / H^ and parameter plane for a phase transition with vacuum energy a — 0.1 (left 
panel) and a = 0.5 (right panel) for Model A. Figure 4 shows the LISA sensitivity regions in the a and parameter 
plane (left panel) and in the 5^ and parameter plane (right panel) for Model A. 
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V. CONCLUSION 



We have analyzed the generation of gravitational waves at the electroweak phase transition, from hydrodynamic 
turbulence in the presence of a helical magnetic field. It is convenient to consider this as a two-stage process. During 
the first stage the effects of helicity are negligible when studying the production of gravitational waves from turbulence. 
The first stage ends when quasi-conservation of magnetic helicity starts to trigger an inverse cascade of the magnetic 
field, that is, a transfer of magnetic energy from small to large scales. The details of the inverse cascade process 
in helical MHD are still under debate and there exist in the literature two different models. In this paper we have 
examined gravitational wave generation in both models of helical MHD turbulence. We have found that, for realistic 
values of the parameters defining the electroweak phase transition, the generated gravitational wave spectrum is 
independent of turbulence modeling (within a factor of order 2). 

The inverse-cascade associated peak frequency always coincides with the Hubble frequency at the moment of 
gravitational wave production. Moreover, as clear from Fig. 1, the main contribution to the total gravitational waves 
energy density is from the second, inverse-cascade stage, even for small values of magnetic helicity. 

Gravitational waves produced via helical turbulence are strongly polarized, since magnetic helicity is a parity-odd 
quantity and is maximal at the end of the first stage [13] ■ LISA should be able to detect such a gravitational wave 



polarization 2^. Also, in contrast to the unmagnetized case, the contribution to the gravitational wave amplitude 



coming from the inverse-cascade stage is large enough at 0.1 mHz to be detectable by LISA. If the electroweak phase 
transition occurs at higher energy (see Fig. 2) the peak frequency is higher, closer to the LISA sensitivity peak, 
which leads to a stronger signal. Our formalism is also applicable to gravitational wave production at an earlier QCD 
phase transition, assuming the presence of colored magnetic fields ^Eu\, or to any other phase transition [6^1 : the peak 
frequency will shift according to changes in Ti, and g^,. 

Finally, we stress that the gravitational wave signal arising from helical MHD turbulent motion exceeds that from 
bubble collisions 0, [13, [IB] and that from hydrodynamical turbulence alone [TO, [lB|- Of course, the strong signal 
estimated here (see also Ref. [ij) assumes initial non-zero (although small) magnetic helicity, so detection of pola rized 
gravitational waves by LISA will indicate parity violation during the phase transition, as proposed in Refs. [20|. 
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